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ABSTRACT 

We present new observational determination of the evolution of the rest-frame 70 and 160 
torn and total infrared (TIR) galaxy luminosity functions (LFs) using 70 \xm data from the 
Spitzer Wide-area Infrared Extragalactic Legacy Survey (SWIRE). The LFs were constructed 
for sources with spectroscopic redshifts only in the XMM-LSS and Lockman Hole fields 
from the SWIRE photometric redshift catalogue. The 70 |xm and TIR LFs were constructed 
in the redshift range < z < 1.2 and the 160 \xm LF was constructed in the redshift range 
< z < 0.5 using a parametric Bayesian and the l/V max methods. We assume in our models, 
that the faint-end power-law index of the LF does not evolve with redshifts. We find the the 
double power-law model is a better representation of the IR LF than the more commonly 
used power-law and Gaussian model. We model the evolution of the FIR LFs as a function of 
redshift where where the characteristic luminosity, L* and the LF normalisation, <f>*, evolve 
as oc (1 + z) Ul and oc (1 + z) a ° respectively. The rest-frame 70 |im LF shows a strong 
luminosity evolution and negligible density evolution out to z = 1.2 with a L = 3.44+g '^g. The 
rest-frame 160 torn LF also showed rapid luminosity evolution with a L — 3.70lg'24 and a D = 

-0.041 o 40 ou t to z — 0.5. The rate of evolution in luminosity and density is consistent with 
values estimated from previous studies using data from IRAS, ISO and Spitzer. The TIR LF 
evolves in luminosity with a L = 3.531q 2 5 and the normalisation evolves as a D — — 0.151q 34 
which is in agreement with previous results from Spitzer 24 \im which find strong luminosity 
evolution and marginal density evolution. By integrating the LF we calculated the co-moving 
IR luminosity density out to z = 1.2, which confirm the rapid evolution in number density 
of LIRGs which contribute ~ 70% to the co-moving star formation rate density at z = 1.2. 
Furthermore, our model suggest that by z = 1.2, ULIRGs are responsible for ~ 10% of the 
IR luminosity density. Our results based on 70 |am data confirms that the bulk of the star 
formation at z — 1 takes place in dust obscured objects. 

Key words: Galaxies: evolution - infrared: galaxies - galaxies: starburst - cosmology: ob- 
servations 



1 INTRODUCTION 

The galaxy luminosity function (LF) and its evolution is an impor- 
tant probe which describes the distribution of galaxies as a function 
of luminosity over the history of the Universe. LFs have been used 
to constrain galaxy formation and evolution models and to quan- 
tify the evolution of the star formation rate (SFR). The infrared 
(IR) LF is essential to understand the amount of energy released by 
reprocessed dust emission from star formation and active galactic 



nuclei (AGN). The discovery of IR lum inous galaxies from ground 
based photometry dRieke & Lowlll972h and by the Infrared Astro- 
nomical Satellite (IRAS) have found them to be locally very rare 
and to only contribute ^5 % to the local IR luminosity density 
JSoifer & Neugebauerlll99lh. The detection of the Cosmic Infrared 
Background (CIB: |Puget et alJI 19961 : iFixsen et al]|l998l) however, 
has shown that roughly half of the energy released in the Universe 
has been absorbed by dust and re-radiated into the IR, which im- 
plies that dust obscured star formation was much more important 
at higher redshifts. 
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Studies using observations performed by IRAS and the ln- 
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frared Space Observatory (ISO) have shown that dusty star forming 
gala xies have undergone st rong evolution, as demonstrated by their 
LFs. lSaunderse t al. ( 1990) construct the 60 um and 40 — 120 um 
far-IR (FIR) LFs based on IRAS observations finding strong lumi- 
nosity evolution (modelled as L*(z) oc (1 + z) aL , where L* is 
the characteristic luminosity and z is the redshift) with a L = 3 ± 1. 
Similar rates of positive evolution (a L = 3 — 5) are seen in LF s 
construc ted from ISO surve ys at 12 um dClements et al. 2001 ). 



90 urn dSeneant et alj | 2004l) and 170 um dTakeuchi etal 



2006) 



IPozzi et alj 1 2004 ) determine the 15 um LF of galaxies from the 



European Large Area ISO Survey (ELAIS) to find that the starburst 
population evolves both in luminosity, with a L = 3.5, and density 
(modelled as <f}*( z ) oc (1 + z) a ° , where 0* is the characteristic 
number density) with a D = 3.8 being consistent with model predic- 
tions of source counts and redshift distribution. 

The sensitivit y and spatial resol ution of the Spitzer Space 
Telescope (Spitzer; IWerner et al. 1 120041) has revolutionised our un- 
derstanding of the evolution of IR galaxies particularly at high 
redshifts (z > 1). Several studies based on Multiband Imager 
for Spitzer (MIPS) 24 um o bservations have been used to con- 
struct the rest-frame 8 um dBabbedge et alj 120061: ICaputi et al.l 
l2007l ; lRodighiero et al.l2010h . 12 urn JPerez-Gonzalez et al l2005l) . 
15um jLeFloc'hetal l 120051). 24 urn teabbedge et all 120061 : 



Ruiopakarn et al. 2010; Rodighiero et al. 2010) and total IR LFs 
dLe Floc'h et alj 20051: [Caputi et al.l 120071 : iMagnelli et alj 120091 ; 
Rodig hiero et aljboid) . These studies have found strong luminos- 
ity evolution with a L ~ 3 — 5 and moderate density evolution out to 
z ~ 1 implying IR galaxies wer e more luminous and numer ous at 
higher redshifts than at z = 0. IPerez-Gonzalez et alj (2005) anal- 
ysed the IR galaxy LF to higher redshifts (z ~ 3) and found that 
the evolution remains constant from z ~ 1.2 to 3. 

All these studies have shown an evolution of IR LFs with look- 
back time and the relative contribution from quiescent (Lir < 10 
L©), luminous IR galaxies (LIRGS: Lir = 10 11 - 10 12 L ) and 
ultraluminous IR galaxies (ULIRGS: Lir = 10 12 - 10 13 L ) to 
the cosmic SFR density. At z < 0.5, the SFR is dominated by 
quiescent galaxies whereas at z ~ 1 LIRGs are responsible for 
~ 50% of the total I R density and dominate the s tar-forming ac- 
tivity beyond z ^ 0.7. Perez-Gonzalez et al. N2005I) show that this 
evolution continues up to z ~ 2.5, and tha t ULIR Gs play a rapidly 
increasing role for z 1.3. ICaputi et alj d2007t) have shown that 
at z ~ 2, around 90% of the IR luminosity density associated with 
star-formation is produced by LIRGs and ULIRGs. Recent results 
from analysis of IR LFs using data from the AKAR1 satellite have 
found good agr eement with the previous studies from I RAS, ISO 
and Spitzer (see lBefhermin et alfeOlOaklGoto et alfcoiol) . 

Most of the previous LF work has been carried out in the mid- 
IR (MIR; A = 8 — 40p.m), while studies of LFs at FIR wave- 
lengths (A = 40 — 200/xm) have been res t ricted to z < 0.3 (se e 
ISaunders et ail Il990l : ISerieant et ail 120041 : iTakeuchi et all l2006t) . 
Studying the redshift evolution of LF at FIR wavelengths is vi- 
tal because the CIB and the SEDs of most IR luminous galaxies 
both peak in this region of the IR spectrum. Furthermore, most of 
the progress in understanding the evolution of IR LFs has utilised 
Spitzer 24 um observations, which are strongly dependent on the 
SED library as at high redshifts the 24 um channel sample shorter 
wavelengths. Further uncertainties in obtaining reliable estimates 
of the bolometric IR luminosities is introduced by redshifting of 
the polycyclic aromatic hydrocarbon (PAH) emission and silicate 
absorption features into the 24 um band at z > 1. 

In this work, we investigate the evolution of the rest-frame 
70um and total IR (TIR) LFs out to z ~ 1.2, and 160 um LF 



out to z ~ 0.5. We use optical and IR data from the XMM- 
LSS and Lockman Hole (LH) regions of the Spitzer Wide-area In - 
fr aRed Extragalact i c surv ey (SWIRE: ILonsdale et al.ll2003ll2004t) . 
In lBabbedgeetalJ 120061) SWIRE optical and IR data were used 
to estimate photometric redshifts (photo- z 's) and construct LFs at 
3.6, 4.5, 5.8, 8 and 24 um over the redshift range < z < 2 for 
galaxies and < z < 4 for optical quasi-stellar objects (QSOs). In 
the present work, the LFs are constructed using spectroscopic red- 
shifts taken from the literature and our spectroscopic follow-up of 
SWIRE selected 70 um sources. Analysis using photo-zs work best 
at high redshifts, where fractional errors of Az/(1 + z) ~ 0.05 
would be acceptable. For example a source at z = 1, the uncer- 
tainty in z would be a tolerable 10%, whereas a source at z — 0.3 
would have an uncertainty of ~ 22%. 

The paper is organised as follows: in Section 2, we present the 
optical and IR data and the spectroscopic redshifts from SWIRE- 
XMMLSS and SWIRE-LH regions used in this study. In Section 3 
we describe the methodology used in calculating the LFs, starting 
with derivations of luminosities, parametrisation of the evolution 
of the LF with redshift and the estimation of SFR with the LF. The 
results of the LFs are presented in Section 4 including estimation of 
the evolution of the integrated SFR and further discussed in Section 
5. We adopt a flat ACDM cosmology model with Ho = 70 km s -1 
Mpc" 1 andf2 A = 0.7. 



2 THE DATA 



2.1 Infrared Data 



The SWIRE survey dLonsdale etal] I2003L |2004|) is one of the 
largest Spitzer legacy programmes covering a 49 deg 2 in six dif- 
ferent fields (ELAIS-N1, ELAIS-N2, ELAIS-S1, CDFS, Lockman 
Hole and XMM-LSS) with both the IRAC (3.6 to 8 um channels) 
and MIPS (24 to 160 um channels) instruments. Typical 5<r sensi- 
tivities are 3.7, 5.3, 48 and 37.7 ^Jy in the IRAC 3.6, 4.5, 5.8 and 
8um bands. For MIPS the 5a limits are 230 /iJy, 20 mjy, and 120 
mjy at 24, 70 and 160 um dSurace et al.l2005l) . The IRAC data were 
processed by the Spitzer Science Ce ntre (SSC) IRAC pipelin e and 
sources extracted using SExtractor dBertin & Arnouts I ll996l) . The 
MIPS 24 um data were processed by the SSC's MOPEX software 
and source extraction was performed using SExtractor. The MIPS 
70 and 160 um data were also processed at SSC using the MOPEX 
software and source extraction carried out through the PRF fitting 
capabilities of MOPEX. Full de tails of the SWIRE data release can 
be found in Sur ace et al. (2005). 

The final data products consist of a cross-matched IRAC and 
MIPS 24 um catalogue and single-band catalogues at 24, 70 and 
160 um. The IRAC and MIPS 24um catalogue consists of sources 
detected with a signal-to-noise ratio (S/N) > 5 in one or more 
IRAC band and their 24 um associations with a S/N > 3. In more 
than 90% of the cases, the sources in the merged catalogue are 
within 1.5" for IRAC channel-pairs and within 3" for MIPS-24 
to IRA C pairs and hence can be considered reliable. ISurace et alj 
d2005l) evaluated the completeness of the IRAC data by comparing 
to deeper data in the ELAIS-N1 region taken as part of the Extra- 
galactic First Look Survey (FLS) programme. The 95% complete- 
ness level was calculated to be at 14, 15, 42 and 56 (ijy in the IRAC 
3.6, 4.5, 5.8 and 8 um bands. The 24 um data is ~ 97% complete 
at 500 p,Jy dBabbedge et"afll2006l) . 

The MIPS 70 and 160um sources single-band catalogues were 
matched to the IRAC and MIPS 24um bandmerged catalogues to 
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produce a SWIRE IRAC and MIPS 7-band catalogue assuming a 
search radius of 9.3" and 19.5" for 70 and 160 \im sources respec- 
tively, which correspond to the telescope limited resolution of the 
MIPS instrument (Vaccari et al., in prep). Thus, almost all of the 
70 and 160 u.m sources are also detected at 24 (jm. The 70 and 160 
\im catalogues are 90% complete down to 14 and 75 mJy for 70 
and 160 \im sources with S/N > 3 (Vaccari et al., in prep). Full 
details on the data processing, completeness and reliability for the 
SWIRE 70 and 160 urn observations will be provided in Vaccari et 
al. (in preparation). 

In this paper we have used data from the latest release of the 
SWIR E photometric redshift catalogue of Ro wan-Robinson et al.l 
(2008). The catalogue contains photometric redshifts for over 1 
million IR sources, estimated by combining the optical and IRAC 
3.6 and 4.5 urn photometry to fit the o bserved SEDs with a 
combination of galaxy and A GN templates (Babbedge et alj|2004l ; 
iRowan-Robinson et alj|2005r) . Our analysis uses multiwavelength 
data from the XMM-LSS and Lockman Hole (LH) regions of the 
SWIRE survey. 



2.2 Optical Data 

Optical photometry is available for > 70% of the SWIRE area 
in at least three of the U, g, r, i and Z ph otometric bands 
jRowan-Robinson et alj|20"08l : iTrichas et alJl2009T) . Spifzer-optical 
cross-identifications (XID) was carried out between the opti- 
cal a nd the IRAC-24ixm catalogues using a searc h radius of 
1.5" teowan-Robinson et alJbOQa . lSurace et alj|2005l) . The cross- 
identification process ensured that each SWIRE source only had 
one optical m atch. Comple t eness and reliability of the XID was in- 
vestigated by ISurace et al.l (2005), which showed that essentially 
the Spfeer-optical XIDs are 100% complete. The requirement that 
the sources be detected at both 3.6 and 4.5 urn at S/N ^ 5 appears 
to eliminate spurious sources effectively and give a high-reli ability 
catalogueue jRowan-Robinson et alj |2005; Sura ce et alj2005h . 

The LH region centred on a= 10 h 45 m , 5 = +57 d 59 m covers 
~ 10.5 deg 2 . Optical photometry is available covering 7.53 deg 2 
in the g, r and z bands were obtained using the MOSAIC cam- 
era on the 4m-Mayall Telescope at Kitt Peak National Observatory 
(KPNO). The 5a limiting magnitudes (Vega) are 25.1, 24.4, an d 
23.7 in the three bands for point-like sources dBertaetal.ll2007l) . 
(/-band photometry is available in a smaller 1.24 deg 2 region to 
5a limiting magnitude (Vega) of 24.1. The data reduction was per- 
formed with the Cambridge Astronomical Survey Unit (CASU) 
pipeline dlrwm&Lewis]|200lh . 

The XMM-LSS field centred on a = 02 h 21 m , S = -04 d 30 m 
and covers 9.1 deg 2 . Optical data is available for 6.97 deg 2 of 
XMM-LSS, which was observed as part of the Canada-France- 
Hawaii Telescope Legacy Survey (CFHTLS) in the u, g, r, i, and z 
bands to magnitude (Vega, 5a for a point like object) limits of: 24.9, 
26.4, 25.5, 24.9 and 23.4 respective^ The photometry was taken 
from IPierre et al. I J20071) . In addition there is 10-band photometry 
(ugrizUB VRI) from the VIMO S VLT Deep Survey (VVDS) pro- 
gramme dLe Fevre et al.ll2004l) covering 0.79 deg 2 and very deep 
5-band photome try (BVRiz) in 1 12 de g 2 of the Subaru XMM Deep 
Survey (SXDS; lFurusawa et alj|2008l) . 



1 See http://www.cfht.hawaii.edu/Science/CFHTLS/ 
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Figure 2. Spectroscopic redshift distribution for the final 70 um sample. 
The 160 u.m sample is a subset of the 70 u.m sample restricted to z < 0.5. 

2.3 Sample Selection 

To study the 70 and 160 [im LFs and their evolution, we select 1.) 
70(im sources with r < 22 mag and S70 > 10 mly; and 2.) 160 
[im sources with r < 22 mag and Sieo > 60 mJy. Furthermore, we 
require the sources to have measured spectroscopic redshifts. The 
adopted 70 and 160 p.m flux limits correspond to 70% complete- 
ness of the SWIRE survey. The total IR LF is derived using the 70 
(jm sample. 

The SWIRE photom etric redshift catalogue of 
IRowan-Robinson et alj d2008l) is the parent catalogue from 
which we select the 70 and 160 urn sources. In LH there are 
4046 and 1276 sources detected at 70 and 160 urn respectively. 
Applying our selection criterion, 2159 70 urn and 584 160|j.m 
sources were selected. Spectroscopic redshifts are available for 
354 70 urn and 166 160 urn sources, which were collated from 
literature. 

In XMM-LSS, 2358 sources have 70 urn detections. We se- 
lected 1606 70u.m sources, of which 299 have spectroscopic red- 
shifts At 160 urn we found 1604 sources with 826 brighter than 
the flux limits stated above. We have spectroscopic redshifts for 
81 160 [im sources. The spectroscopic redshifts in XMM-LSS 
were obtained from our spectroscopic follow-up programme con- 
ducted over 6 nights between November 2008-November 2009 at 
the WHT using the AF2/WYFFOS inst rument. The detai ls of the 
spectroscopic follow-up are presented in lPatel et al.l J201 lh . 

We choose to construct the 70 (im and the total IR LFs in the 
redshift range < z ^ 1.2 and the 160 [im LF in the redshift range 
< z ^ 0.5. Thus, the final total sample of galaxies considered 
for the 70 [im and total IR LFs consists of 634 (343 in LH and 
291 in XMM-LSS) objects and for the 160 LF consists of 221 
(157 in LH and 64 in XMM-LSS) objects. We show in Figure[T]the 
r-band magnitude as a function of the 70 u.m flux for all sources 
in the LH and XMM-LSS regions (filled blue circles) and the final 
spectroscopic sample (filled red circles). 

A summary of the parent catalogue and the final sample is pre- 
sented in TableQ] We display the spectroscopic redshift distribution 
of the final 70 [im sample in Figure[2] while the final 160 urn sam- 
ple has the same redshift distribution for z ^ 0.5, since all the 160 
[im detected sources are also detected at 70 (im. The redshift dis- 
tributions of the two fields is different because of large scale struc- 
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Figure 1. Left: r-band magnitude against 70 urn flux. Blue dots are all 70 urn detected sources with an optical counterpart in the LH and XMM-LSS. Filled 
red circles are the sources in the final 70 |j.m sample. The vertical line marks the 70 |j.m selection limit of 10 mjy and the horizontal line is at r = 22. Right: 
r-band magnitude against 160 |xm flux. Blue dots are all 160 |^m detected sources with an optical counterpart in the LH and XMM-LSS . Filled red circles are 
the sources in the final 160 |j.m sample. The vertical line marks the 160 u.m selection limit of 60 mJy and the horizontal line is at r = 22. 



Table 1. Summary of the sample selection used to study the 70 and 160 \im LFs. 



Field 


a (J2000) 


8 (J2000) 


Survey area 




70 urn sources 






160 u.m sources 










(deg 2 ) 








iV sp d 


N iet a 


AV 


N sp c 


















z < 1.2 








z < 0.5 


LH 


10 h 45 ra 22 s 


+57 d 59 ra 05 s 


7.53 


4046 


2159 


354 


343 


1276 


584 


166 


157 


XMM-LSS 


02 h 21 ra 20 s 


-04 d 30 m 00 s 


6.97 


2358 


1606 


299 


291 


1604 


739 


81 


64 



"Number of sources detected; ^Number of sources selected using the selection criteria; c Number of sources with spectroscopic redshift; d Number of sources 
with spectroscopic redshift below 1.2; e Number of sources with spectroscopic redshift below 0.5 



tures in the XM M-LSS field where the cluster distribution peaks 
around z = 0.3 l lPacaud et al.l2007ri . which also corresponds to the 
peak in the redshift distribution of our sample. 



3 INFRARED LUMINOSITIES 

To determine the 70, 160 \im and total IR LFs from our sample, we 
need to derive the rest-frame 70 and 160 \im and the total IR lu- 
minosities. In order to do this, we model the SEDs for each source 
follow ing the method described in iRowan-Robinson et alj ( I2005L 
2008) used for the SWIRE photometric redshift catalogueue, using 
optical (at least 3 of the 5 optical U, g, r, i and Z bands) and IR pho- 
tometry (Spitzer IRAC 3.6 - 8 \im and MIPS 24 - 160^m bands). 
The SED fitting follows a two-stage approach, by first fitting the 
optical to near-IR (U to 4.5 itm) SED using the s i x gala xy and three 
AGN templates used bv lRowan-Robinson et alj 1(2008). 

We calculate the IR excess by subtracting the galaxy model 
fit from the 4.5 to 24|im data. We then fit the IR excess, 70 and 
16 0yim (for 83 70^tm sources) data points w ith the IR template 
of IRowan-Robinson et al] d2004l . |2005| . 120081) . The IR templates 
are derived from radiative transfer models dependent on interstel- 
lar dust grains, the geometry and the density distribution of dust. 
The IR templates are: 1.) IR 'cirrus': optically thin emission from 
interstellar dust illuminated by the interstellar radiation field; 2.) an 
M82 starburst; 3.) a more extreme Arp220-like starburst and 4.) an 
AGN dust torus. We also allow the sources to be fit by a mixture of: 
1.) M82 starburst and cirrus, 2.) M82 starburst and AGN dust torus 



14F^~^ — — 1 — " — i -1 — 1 — ■— i — ■— ' — ' — — ■ — i — ■ — ■— i 

: HLIRGs : 




: * 1 ''at. : 

10 T □ K Arp220 - 

: + + M82 : 

Cirrus I 

AGN Dust Torus I 

I Q,; M82 + CIrrus mixture I 

9 T'4 □ M82+AGN Dust Torus mixture ~ 

•Arp220+AGN Dust Torus mixture : 

2L, , , i , , , i , , , i , , , i , , , i , , , [ 

0.0 0.2 0.4 0.6 0.8 1.0 1.2 

Z sp 

Figure 3. Total IR luminosity (Lir ) as a function of redshift for 70 u.m 
sources with flux above 9 mJy. The horizontal lines mark the divisions 
for LIRGs (10 11 Lq < L IR < 1O 12 L ), ULIRGs (10 12 L Q < L IR < 
10 13 L Q ) and HLIRGs (L IR > 10 13 L Q ). 



and 3 .) Arp220 and AGN dust torus to properly represent the IR ex- 
cess jRowan-Robinson & Crawford! fl989l :lR owan-Robinson et al. 
l2005h . 

From the SED fitting we find all but 12 of the z < 1.2 sample 
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Figure 4. Left: Rest-frame 70 [im luminosity (L70) as a function of redshift. The lines indicate as a function of redshift the 70 [im rest-frame luminosity 
corresponding to an observed 70 u.m flux of 10 mJy with the 4 IR templates. Right: Rest-frame 160 |rm luminosity (L160) as a function of redshift. The lines 
indicate as a function of redshift the 160 u.m rest-frame luminosity corresponding to an observed 160 u.m flux of 60 mJy with the 4 IR templates. 



are fit with a galaxy template in the optical, which implies that the 
contamination by AGNs is low in our sample. Thus, we choose to 
include all sources in our analysis and do not separate the galaxies 
from AGNs. We illustrate in Figure [3] the distribution of total IR 
luminosity (estimated by integrating the SED models between 8 - 
1000 (xm) as a function of redshift for the z < 1.2 sample, with 
symbols representing the best fitting IR SED. The figure shows 
that in the IR the entire sample is dominated either by the starburst 
(Arp220 or M82-like) or the cirrus component, which is consistent 
with previous stu dies dTrichas et al. 2009; Svmeonidis et alfeoid : 
IPatel et al.l201lh . Figure[3]also shows that most of the 70 pm sam- 
ple are LIRGs and none of our sources have a pure AGN dust torus 
dominated IR SED. 

We show in Figure [4] the rest-frame 70 [im luminosity (L70) 
distribution as a function of redshift (left panel) and the rest-frame 
160 (xm luminosity (L160) distribution as a function of redshift 
(right panel). The 70 and 160 [im luminosities were computed us- 
ing K-corrections derived from the best fitting model templates. 
The uncert ainty in the 70 and 160 [im and Lir are typically 0. 1 dex 
as stated in IPatel et al.l f2oTTh . The SED model templates are also 
used to determine 1/V max , which is the volume corresponding to 
the maximum redshift at which a source could be detected by the 
survey given the optical and IR limits set in Section l2~3l 



4 METHODOLOGY 

In this section, we describe the methods used to determine the red- 
shift evolution of the rest-frame 70, 160 pm and TIR LFs, derived 
using the spectroscopic sample described in Section 12.31 and us- 
ing the optical and IR SED model of each source described in 
the previous section. The LF, denoted as <3?(-L) dlogL, defines 
the number of objects per comoving volume within a luminosity 
range logL, logL + dlog L. We construct the 70 (xm and total IR 
LFs in the redshift bins [0-0.2], [0.2-0.4], [0.4-0.8] and [0.8-1.2] 
and the 160 pin LF in the redshift bins [0-0.2] and [0.2-0.5 ] us- 
ing the l/y m ax method dSchmidJ[l968llAvni & Bahcaiilll980l) . We 
also develop a Bayesian method to study the evolution of the 70 
and 160 [im and the TIR LFs. In our LF estimation methods we 



accurately model the different completeness and selection effects 
affecting our data. The selection and incompleteness functions are 
described in Section l4~T1 and the binned and parametric LF methods 
are discussed in Section |4~2l and Section |4~3"l respectivelv. 



4.1 Selection and Incompleteness Function 

The data set we have used to estimate the FIR LFs was selected by 
imposing multivariate flux limits to determine whether a source is 
included in our sample. In addition we require each of the sources 
to then have a spectroscopic redshift. Here we use the 70 pm LF 
selection functions as an example but the general procedure is the 
same for the total IR and 160 (im LFs. 

The 70 pm LF is calculated by selecting sources in the SWIRE 
catalogue that have S70 > 10 mJy, r < 22 mag and measured 
spectroscopic redshift. Therefore our first selection function con- 
siders the probability that a source of a given 70 [im flux (which 
is a function of luminosity, L, and redshift z) is detected by the 
survey. This is denoted as p(det|L, z), and has been determined by 
Vaccari et al. (in preparation) for each field in the SWIRE survey 
using Monte Carlo simulations. The second selection function con- 
siders the probability that a given 70 [im source with luminosity, L 
at redshift z is associated with an optical counterpart at r < 22 mag 
and is defined as p(r < 22\L, z). This was quantified by taking the 
SWIRE catalogue and constructing a source count distribution for 
all sources detected at 70 pm and then constructing a similar source 
count distribution for all 70 (im sources with r < 22 mag and the 
ratio is used to estimate p(r < 22\L, z). 

Finally, we take into account the probability that a given 
source to be characterised with a spectroscopic redshift. The spec- 
troscopic incompleteness is primarily determined by the r band 
magnitude and therefore they are added as weights in the com- 
putation of the LF (see Section l4~2l and Section l4~3l . We define 
the weights, w, as the inverse of the spectroscopic completeness. 
The spectroscopic completeness is determined by constructing the 
source distribution for all sources with a measured spectroscopic 
redshift as a function of r magnitude and dividing this by the source 
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count distribution for all 70 (J.m sources with S70 > 10 mJy and 
r < 22 as a function of r magnitude. 

The selection functions were estimated for LH and XMM-LSS 
regions separately and used in constructing the 70, 160 pm and 
total IR LFs. The selection function, defined as p(selected|L, 2), 
are combinations of the two selection functions: 



p(selected|L, 2) = p(det|L, z)p{r < 22\L, z) 



(1) 



The selection function and the spectroscopic incompleteness were 
used to modify the LF methods discussed in the following section. 



4.2 Binned Estimates 

We use the 1/Kn ax method as the binned estimate which has the 
advantage that it allows direct computation of the LF from the 
data, without any parametric dependence or model assumption. We 
divided the sample into redshift bins selected to ensure adequate 
numbers of galaxies in each bin. For each redshift bin, the LF is 
given by: 



$(L)dlogL= \ Y,m X T7~ J 

\ . vmax,i J 



(2) 



where V max ,i is the comoving volume out to which the rth galaxy 
could be observed, Wi is the inverse of the spectroscopic incom- 
pleteness of the ith galaxy. The comoving volume, V maXj , is: 



dV 

nfselectedlLi, z) — — dz, 
dz 



(3) 



where dV/dz is t he differential comoving element per unit solid 
angle (Hogg 1999). Here 2 max .; corresponds to the maximum red- 
shift at which the source could be detected by the survey given the 
optical flux limit (r < 22) or the IR flux limits (S70 > 10 mJy 
or Sieo > 60 mJy) and 2 m i n is the lower limit of the redshift bin. 
Zmax,i was determined by using the optical/NIR and IR SED model 
of each source. The volume element integral in the LF calculation 
was weighted by the selection function to correct for the selection 
biases that are inherent in the spectroscopic catalogueue. The IR 
SED models described in Section [3] are used in computing the se- 
lection function. 

The associated rms error is given by: 



E'< 



V 2 



(4) 



In order to accurately determine the uncertainty in the LF error 
we use Monte Carlo bootstrapping analysis to randomly resample 
the final spectroscopic catalogueue to generate 1000 realisations; 
each of these is analysed as described above and the root mean 
square (RMS) of the results is quoted as the error. 



4.3 Parametric Bayesian Method 

We use a Bayesian approach to determine the parametric LF, which 
requires prior knowledge of the appropriate functional form of 
the LF. As is the c a se for maximum lik elihood (ML) methods 
dSandageet alii 19791 : iMarshall et all l983) the parametric LF has 
the advantage over the \IV ma x. method in that it is insensitive to 



any local clustering effect whereas the 1/Vm ax LF assumes a uni- 
form number density throughout the observed volume and there- 
fore is vulnerable to density inhomo geneities present in the sur- 
vey JWang & Rowan-Robinsor] l2010). The advantage of using a 
Bayesian method to estimate the LF parameters over ML methods 
is that ML methods do not provide an estimate of the LF normalisa- 
tion, which is often chosen to make the expected number of sources 
detected in a survey equal to the actual number of sourced detected. 
In addition, the confidence intervals on the LF parameters are de- 
rived assuming they have a Gaussian distribution whic h is not nec- 
essari ly a good approximation for small sample sizes (K elly et all 



In order to carry out a Bayesian analysis we first need to de- 
fine the likelihood function, p({d}\{8}), which is the probability 
of observing the data, {d}, for a given LF model, that is described 
by some parameters {8}. To do this we first define the probabil- 
ity of finding a source of a specific luminosity Li in the range 
log L, logL + d log L at a redshift Zi in the range z, z + dz as: 



p(L,z\{8}) = 



$(L, 2|{6>})p(selected|L, 2) dV 



A d2 ' 

where A is the expected number of sources and is determined by: 



(5) 



E 



dV 

$(L,2|{0})p(selected|L,2)dlogL-jjd2 (6) 



The sum is taken over the fields present in our survey and the inte- 
grals are taken over all possible values of redshifts and luminosities. 

We now write the likelihood function as the probability of ob- 
serving N objects, each with Li and Zi, drawn from the model LF 
as: 



P ({d}m)=p(N,{Li,Zi}\{e}) 

= p(N\{e})p({L i ,z i }\{0}) 



(7) 



where p(N\{9}) is the probability of observing TV objects given 
the model LF and p({Li, Zi}\{9}) is the likelihood of observing a 
set of Li and z% given the model LF. We assume that the number of 
sources detected follows a Poisson distribution, where the expected 
number of detectable sources, A, is given by Equation[6] Thus the 
likelihood function is written as: 



p(N, {Li, Zi}\{8}) 



■JJp(Li,2,|{fl}) 



n 



A"e- 

7V"! 

i— 1 

A e 

Hvr x 

$(L, 2|{0})p(selected|L, 2) dV 



A 



dz 



(8) 



(9) 



We further modify the likelihood function by including the 
spectroscopic incompletenes s by introducing a weighting factor , 
w/(w), for each object (see IZucca et al.lll994l : lllbert et al]|2005t 
lAird et al.ll2008h . The weights, added as exponents of the indi- 
vidual source likelihoods p({Li, Zi}\{9}) artificially reduce the 
size of error estimates. Therefore they are balanced by the aver- 
age weight ((w)) which do not effect the best fitting parameters 
dAird et all 120081) . Furthermore, N, the total number of objects 
is now J^lOj, which gives the effective number of sources cor- 
rected for the spectroscopic incompleteness. Therefore the likeli- 
hood function in Equation|9]is: 
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p(N,{U, Zl }\{e})(xX^^e- x x 

if 

n 



$(L, z|{6>})p(selected|L, g) dV^ \ ~ 
A dz J 

(10) 



Note that in the absence of the spectroscopic incompleteness 
weights, Wj, Equati on [Tol reduces to the form presented in 
lMarshalletal](ll983h . 

We perform Bayesian inference by combining the LF with 
a prior probability distribution, p({9}) to compute the posterior 
probability distribution, p({6}\{d}), given by Bayes' theorem: 



P({0}\{d}) = 



P({d}\{0})p({0}) 
Jp({d}\{6})p({8})d6> 



(11) 



where the denominator is the Bayesian evidence and is determined 
by integrating the likelihood over the prior parameter space. 

For parameter inference, the Bayesian evidence serves to nor- 
malise the posterior distribution and is vital for Bayesian model 
comparison. Calculating the Bayesian evidence is computationally 
expensive since it involves integration over n-dimensions for an 
n parameter LF model. Therefore we use standard Markov chain 
Monte Carlo (MCMC) methods to perform a random walk through 
the parameter space to obtain random samples from the poste- 
rior distribution. We employed the Met ropolis-Hastings algorithm 
(Metr opolis et al.|[l953l : iHasting j| 1970l) . in which a proposal dis- 
tribution is used to guide the variation of the parameters. The al- 
gorithm uses a proposal distribution which depends on the cur- 
rent state to generate a new proposal sample. We accept a step 
if the probability of the model given the new parameter values is 
higher and also at random intervals when the probability is lower 
in order to allow the fit to proceed downhill to avoid local minima 
dPtak et alj2007l : lKellv et all2008h . 

We assume a flat prior distribution for each parameter and ig- 
noring the normalising factor, the Bayesian evidence, the posterior 
distribution for each parameter is then given by: 



P ({0}|{d})<xp({d}|{0}) 



(12) 



We produce three chains for each analysis of at least 2 x 10 iter- 
ations and adjust the parameter step sizes to achieve an acceptance 
ratio in the rang e 0.3 - 0.5. Finally w e calculated the convergence 
R statistic from lGelman et al.l ( l2004) . which should be sC 1.2 if the 
chain has converged. For all the parameters the R value was < 1.1. 



4.3.1 Luminosity Function Models 

The first IR LF constructed from IRAS observations showed an 
excess in the num ber of galaxies at the high luminosity end 
dSoiferet al.lll987h from the value expected from the Schechter 
function jSchechtej[l976l) . ISoifer et al.l Jl987b fit a double power- 
law model to the IRAS 60 (am LF which is adopted in our work. We 
use a continuous double power-law model given by: 



*{L\{8}) = , 



+ 



(13) 



Alternatively ISaunders et~aLl dl990h fit the IRAS 60 LF with a 
combination of a power-law and log-normal LF model: 



*(L\{e})=4>* — 



(l-c) 



exp 



1 1 2 

"2^ l0g 



(14) 



In both models, L* is the characteristic luminosity, <j>* is the 
LF normalisation and a is the power-law of the faint-end of the LF. 
P in Equation[T3]is the power-law index of the bright-end of the LF 
and a in Equation [l4]gives the range over which the LF drops off. 

We estimate the parameters of the double power-law model, 
defined model 1, and the power-law and log-normal model, de- 
fined model 2. Uniform prior probability for all parameters were 
assumed with limits as follows: — l^a^2, O^ff^l, 
1 sC P sC 6, 8 < logL* sC 12 and -1 < log0* sC -4. 



4.3.2 Luminosity Function Evolution Model 
Sever al studies (se e Le Floc'h et al.l 20051; Perez-Gonzalez et al 



| 2005 |: ICaputi et alj 120071 : iMagnelli et al.l l2009l : iRuiopakarn et al 
120101) of the IR luminosity functions, mainly at 24 [im have shown 
that out to z ~ 1.2, the while the luminosity and the number den- 
sity of the LF evolves, the shape remains the same, that is a, P (in 
Equationll3t or a (in Equationll4t do not change. We parameterize 
the evolution of the FIR LFs in luminosity only: 



L 

W) 



iw 



(15) 



where f(z) — (1 + z) aL . Therefore in addition to the parameters 
listed in the previous section, a L , the luminosity evolution power- 
law index was also constrained. Uniform prior probability distribu- 
tions were assumed for q l with limits: 1 ^ a L ^ 7. Thus, for each 
LF model (double power-law and power-law and log-normal) five 
parameters were estimated. 



5 RESULTS 

In this section, we present our determinations of the FIR LF and 
compare them to previous studies. We construct the rest-frame 70 
\im and the TIR LFs in redshift range < z ^ 1.2 and the rest- 
frame 160 |J.m LF in the redshift range < z ^ 0.5. Finally using 
the TIR LF we derive the co-moving IR luminosity density and the 
cosmic star formation rate density as a function of redshift in the 
range < z ^ 1.2. 



5.1 Evolution of the Rest-frame 70 fim Luminosity Function 

In Figure [5] we show the rest-frame 70 (im LF constructed using 
the 1/Vmax (black filled circles) and the parametric Bayesian (solid 
red and blue lines) methods described in Section |4] The paramet- 
ric LFs are displayed using the best-fit (posterior mode) parameters 
given in Table [3] We also display the z = LF (dashed red and 
blue lines), which shows the rapid evolution of the 70 vim LF, when 
compared with the binned and evolved parametric LFs. The para- 
metric 70 |J.m LF uses k-correction according to the M82 starburst 
SED. 

Examination of the binned and parametric LFs in Figure [5] 
shows excellent agreement across all redshift bins. Comparison 
of the two LF models indicates that they are almost identical ex- 
cept at the brightest luminosities (L70 > 10 12 L©), where the 
power-law and log-normal LF model has a rapid drop off when 
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log L„ (L s ) log L„ (L e ) 



Figure 5. 70 p LF split into redshift bins < z < 0.2, 0.2 < z sg 0.4, 0.4 < z < 0.8 and 0.8 < z < 1.2. The 1/V max LF is shown as filled black circles, 
the solid blue line is the double power-law LF and the solid red line is the power-law and log-normal LF. The parametric LFs are displayed using the best-fit 
parameters listed in Tableland evaluated at z = 0.1, 0.3, 0.6 and 1. The dashed red and blue lines are the 2 = LFs. 



log L 70 (L Q ) < z < 0.2 



log* (Mpc -3 ) 
0.2 < z < 0.4 0.4 < z < 0.8 



0.8 < z < 1.2 



8.4 
8.8 


-1 Q6 +0M 
— 1.20 

1 OQ+016 




9.2 


9 oq+0.15 
z.oo_ 22 




9.6 


9 oq+0.08 




10.0 


_o KQ+0-06 
z - oy -0.07 




10.4 


9 orr + 0.05 

z -°'-0.06 


-2 5Q+ - 20 
z.oa_ 40 


10.8 


-3.491°;?? 


_ 2 QQ+ 07 
z - 9y -0.08 


11.2 


-4.7il°;^ 


q cc+0.07 


11.6 


r 9S+0-31 
,J - zo -5.28 


_ 4 40+0.13 
^■^ 9 -0.19 


12.0 




d - oz -5.62 



12.4 
12.8 
13.2 



-3.25 
-3.80 
-4.60 
-5.14 



+0.12 

0.16 

+0.10 

0.14 
+0.13 
0.18 
+ 0.16 
0.25 



A 40+O.85 

0.20 
0.37 
+0.25 
0.63 
0.30 
6.36 



-4.52 
-4.70 
-6.36 



Table 2. 1/V m ax 70 |^m luminosity function values. 





70 urn LF 


160 \im LF 


Total IR LF 


Parameter 


(a) (b) 


(a) (b) 


(a) (b) 



logL*(L ) 10.35^ 9.53i°;i° 10.21+"° 9.85t°; 4 * 10.6li°;« 9.71™ 



,+0.10 j 4n +0.10 4n +0.17 1 n7 +0.14 „ o.+O.lO , oa+0.09 

"•"'-0.13 1# * u -0.09 u -* u -0.21 1 ' u '-0.58 U ' J ^ 

2.00+nl? - 3.79+nll - 2.09 



+0.13 o 7Q+0.36 9 no+0.14 

0.12 "-' a -0.23 z - uy -0.10 

+0.07 n 9 Q +0.08 n 00+O.O6 



0.66^;- - 0.29lg;S - 0.63l . 04 

-2.68t°;°° -2.32±°;S£ -2.841°;}! -2.43±g:iS -2.e8±g;ig -2.26i^ 2 

q OQ+0.12 q 41+0.18 J- 7 q+0.30 C £, + 0.28 q O9 + 0.24 q 89 + 0.25 

o.oa_g 2 2 q 25 " J -' o -0.62 J -' JO -0.23 °- oz -0.20 °- oz -0.16 



Table 3. Best-fitting parameters for of the 70 and 160 |j.m and TIR LFs determined using the parametric Bayesian method. The errors include 68% of the 
posterior probability, (a) double power-law LF model and (b) power-law and log-Gaussian LF model. 
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LF evolves rapidly in luminosity (a L > 3). Studying the LF evolu- 
tion of each SED component will be further investigated in future 
studies (see Section [6~3l >. 



5.2 



Evolution of the Rest-frame 160 /xm Luminosity 
Function 



Figure 6. The LLF (0 < z < 0.2) at 70 urn determined Using 1/Vmax 
(filled black circles) and the parametric Bayesian method (solid and dashed 
black lines) evaluated at z = 0.1. The solid black line is the double- 
power law model and the dashed black l ine is the power - law an d log-normal 
model. Binned IRAS 60 |xm LF from ISaunders et all Jl990h is shown in 
filled blue squares and the parametric 60 p.m LF from Takeuc hi et al 
(solid red line) and Wang & Rowan-Robinson ( 20 uj) (solid green line) are 
corrected for the bandpass differences assuming an M82 type starburst SED. 



The rest-frame 160 |J.m LF constructed using the 1/Knax (filled 
black circles) and the parametric Bayesian (solid red and blue lines) 
methods is shown in Figure [7] The parametric LFs are displayed 
using the best-fit parameters given in Table [3] and the 1/Vmax val- 
ues are listed in Table [4] The z = LF is also displayed (dashed 
red and blue lines), which show the rapid evolution of the 160 \im 
LF, when compared with the binned and the evolved parametric 
LFs. The parametric Bayesian 160 \im LF was calculated using a 
k-correction given by a mixture of M82 starburst (50%) and cirrus 
template (50%). 

The binned and parametric LFs shown in Figure |7]are consis- 
tent with each other except at Li6o > 10 L where both the para- 
metric LFs show a faster drop-off. In fact the double power law LF 
model shows a much better agreement with the l/V m &x LF in both 
redshift bins at Li6o > 10 11 Lq implying that the double power 
law LF model may be a better description of FIR LFs. The 160 (xm 
LF evolves in luminosity with a L = 5.73ig'g2 f° r model 1 and a L 
= 5.531q 23 f° r model 2, which is stronger than the value found 
for the evolution of the 70 [ i.m LF or from Spitzer 24 \im stud- 
ies (seelLe Floc'h et al . 2005; Babbed ge et al.l200 6; Magnel li et al.l 
r2009l : lRuiopakarn et alj20ld) . lTakeuchi et alj ( 120061) however find 

for model 1 and; a = 1.40± o ;$ for model 2 as observed in"fhe that the IS0 170 V™ LF evolves with «>- = 5 -°± 2 oi> being entirely 

consis tent with the evolution of 160 |i.m LF. Recently, IDve et alj 
d2010h used data from the H-ATLAS survey to find that the rest- 
frame 250 p.m luminosity density evolves at a rate proportional to 

7 1 + 21 

(1 + z) ' - 1 4 to z ~ 0.2 in agreement with the evolution of 160 
and 170 p.m LFs. 

Our results show that the 160 [im LF evolves more rapidly 
than the 70 (im LF suggesting that 'cooler' galaxies evolve more 
rapidly than 'warmer' galaxies (the 160 |i.m LF was determined 
using the k-correction given by a mixture of M82 starburst and 
cirrus template whereas the 70 \im LF was determined using the 
k-correction given b y an M82 starburst temp late). This is in keep- 
ing with the work of Svmeonidis et al. I J201 l|) who find that 'cold' 
galaxies evolve more rapidly than 'w armer' galaxies over a period 
of 0.1 < z < l. lDunne et al.l J20TH) have shown that the evolution 
of the 250 |i.m LF out to z = 0.5 is driven in part by evolution in 
the dust mass and an increase in the luminosity or space density of 
cooler galaxies. They conclude that the evolution of the dust mass 
points to an enhanced supply of gas for star formation at earlier 
cosmic epochs. 



compared with the double power-law model. However, since there 
are no data points at these brightest luminosities, both LF mod- 
els are a good description of the evolution of the 70 (im LF. The 
shape of the faint-end of the LF is also consistent; a = 0.371q^3 

-0.09 tr ' ' ' ^ '" k 

< z ^ 0.2 redshift bin. The luminosity evolution parameter, 
a L , which shows that the characteristic luminosity evolves rapidly 
as a function of redshift is almost identical for the two LF mod- 
3 to.22 and »l = 3.41±°; 2 * for model 1 and 2 re- 



3.39Z 



els (q l = 

spectively). The best fit evolutionary parameter is consistent with 
studies performed at several IR wavelengths, in particular at 24 
yim, which find strong luminosity evolution with q L = 3 — 5 (see 
Le Floc'h et alj|2005l; iBabbedge et al.ll2006r . iMagnelli et aT]|2009l ; 



Ruiop akarn et al. 20 id) . We find that both parametric models can 



describe the evolution of the 70 \im LF. 

We shown in Figure [6] the SWIRE 70 \im local luminosity 
function (LLF) constructed in the redshift bin < 2 ^ 0.2, 
is compared to the bandpass corrected (a ssuming an M82 star- 
burst SED) IRAS 60 yun LFs determined by ISaunders et all dl990h 
(filled blue square), ITakeuchi et all J2003b (solid red line) and 
Wang & Rowan-Robinsonl ( feOlOh Tsolid green line). Comparison 
of the LFs shows that the overall shape of the parametric and 
binne d 70 ytm LF is well m atched to the bandpass corrected 60 (xm 
LFs. ITakeuchi et ail d2003l) IRAS 60 yim LF report, log L* (L ) 
= 9.07 ± 0.09, a = 1.23 ± 0. 04, a = 0.72 ± 0.01 and log** 
(Mpc -3 ) = -2.05 ± 0.05 while lWang & Rowan-Robinsotj feoicb 
determine, logL* (L ) = 9.10, a = 1.29, a = 0.72 and log</>* 
(Mpc -3 ) = —2.05. At the faint-end, the 70 |i.m LF is steeper with 



We co mpare in Figure 



\xm LF from Takeuchi et al 



-v+0.10 



5+0.07 



agreement with the value determined from IRAS surveys. 

We have used the k-correction given by the M82 starburst SED 
template in order to determine p(selected|L, z) in the likelihood 
function because a large fraction of the 70 \im population are fit 
with this template (see Section[3]and Figure[3}. Although using a 
k-correction of a pure Arp220 starburst or a pure cirrus SED tem- 
plate results in different values for a L , we still find that the 70 (im 



8] our 160 vim LLF with the ISO 170 
(2006) constructed using 55 galaxies 
at z < 0.3. The 170 \im binned LF shows a reasonable agree- 
ment within the error bars with the 160 |i.m LF except at L160 
< 10 10 L where there is a significant discrepancy between the 
two LFs, which is m ost likely due to the sm all number of objects in 
the 170 (i.m sample. ITakeuchi et all (2006) use the power-law and 
log-normal LF model and conclude that the parametric form (solid 
red line in Figure [8} underestimates the bright end of the LF, sim- 
ilar to the findings reported here. The combination of these results 
indicates that the double power-law model is the favoured analytic 
form for FIR LFs. 
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Figure 7. 160 |^m LF split into redshift bins < z ^ 0.2 and 0.2 < z ^ 0.5. The 1/Vmax LF is shown as filled black circles, the solid blue line is the double 
power-law LF and the solid red line is the power-law and log-normal LF. The parametric LFs are displayed using the best-fit parameters listed in Tableland 
evaluated at Z. = 0.1 and 0.5. The dashed red and blue lines are the z = LFs. 
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Table 4. The 1/V max 160 Aim LF values. 



5.3 Total Infrared Luminosity Function 

In this section, we compute the total IR luminosity function of 
galaxies using the observed 70 data. The total IR luminos- 
ity of a star forming galaxy provides a direct estimate of current 
star formation activity because the IR emission is the reprocessed 
UV/optical radiation produced by young stars. Lir can be con - 
verted to SFR using the relationship provided by Ken nicutJ dl998l) : 

SFR (M yr _1 ) = 1.72 x 10" 10 Lir L , (16) 

where Lir is estimated by integrating for each source the best fit 
SED in the interval 8 — lOOO^m (see Section[5](. We can then use 
Equation [16] and the total IR LF to derive an estimate of the IR 
comoving energy density and the cosmic star formation rate den- 
sity (CSFRD) up to z ~ 1.2 and compare these results with other 
CSFRD calibrators. 

Previous studies of the the total IR LFs have largely been 
conducted at MIR wavelengths, with the greatest progress made 
at the Spitzer 24 \im band where MIPS is most sensitive (see. 
iLe Floc'h etai]|2005l ; IPerez-Gonzalez et al.ll20ol iBabbedge etaD 



2006j : ICaputi et alj|2007l : iBethermin etaill2010bl : iRodighiero etaU 
20101) . The work presented here follows a similar approach to 
others in the literature, which relies on the conversion of L„ to 
Lir , usually calibrated u sing SED templates of l ocal IR galax- 
ies l Chary & Elba z 2001) or semi-empirical SE Ds dDale & Helou 
20021; lLagache et al.l 12003). Several auth ors lElbaz et al.l 120021 ; 
Appleton e t alj |2004|) and in paritucular iBavouzet et al.l ( 120081) . 
have shown that the extrapolation of the local L„ — Lir remains 
reliable up to z = 1.1 for LIRGs and 2 ~ 2 for ULIRGs. Further- 
more, the wavelength closest to the pea k of far-IR emission p ro- 
vides the most accurate estimator of Lir dBavouzet et alj|2008 



We us e the SED templates of iRowan-Robinson et alj d2004l . 
120051120081) and therefore to check the consistency of ourL7p — Lir 
correla tion, we compared our result s with those of IBavouzet et al.l 
(2008) and lSvmeonidis et all J2008I) . In Figure |9]we show the re- 
lationship between log L70 and log Lir for all sou rces at z < 1.2 
in our sample (black solid lin e) with the relation s of IBavouzet et al.1 
(2008) (green solid line) and lSvmeonidis et al.l d2008) (blue solid 
line). The three correlations show excellent agreement with a mean 
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Figure 8. The LLF (0 < z < 0.2) 160 p LF determined using 1/V ma x 
(rilled black circles) and the parametric Bayesian method (solid and dashed 
black lines) evaluated at z = 0.1. The solid black line represents the 
double-power law model and the dashed black line is the power-law and 
Gaussian model. The parametric (red solid line) and non-p arametric (red 
filled squares) ISO 170 urn LFs from Tak euchi et al.l ^2006) are shown for 
comparison. 
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Figure 9. Correlation between rest frame 70 urn luminosity (log L70) and 
total IR luminosity (log Lir). The black solid line is our linear relation- 
ship determined using a Bayesian line ar regression method. W e also show 
for comp arison the linear relation from Symeonidis et al. (2008) (blue solid 
line) and Bavouzet et al. 12008) (green solid line). 



scatter of ~ 0.1 dex between our correlation and the other two in 
the luminosity range 9.5 < logL7o < 12.5. Therefore we choose 
to use our L70 — Lir correlation in constructing the total IR LF. 

To calculate the 1/V max LF, we follow the same method as 
described in Section l4~2l and use the 70 (im selection function for 
each source. For the Bayesian analysis we use our linear relation- 
ship between L70 — Lir to convert Lir to L70 to determine the se- 
lection function p(selected|L, z) in the likelihood function. Thus, 
the Bayesian analysis is dependent on the values of the linear re- 
lationship shown in Figure [9] We perform a test by estimating 
the parameters of the local total IR LF using the correlations of 



iBavouzet etal] J2008h and lSvmeonidis et alj d2008l) and found all 
parameters to be almost identical. 

5.3.1 Evolution of the Total Infrared Luminosity Function 

In Figure [TO] we display the TIR LF constructed using the 1/Knax 
(filled black circles) and the parametric Bayesian (solid black line) 
methods. The parametric LFs are displayed using the best-fit pa- 
rameters given in Tableland the 1/V ma x values are listed in Table 
[5] The z — LF is also displayed (dashed black line), which show 
the rapid evolution of the TIR LF, when compared with the binned 
and the evolved parametric LFs. The parametric TIR LF was cal- 
culated using the k-correction given by the M82 starburst SED. 

Comparison of the 1/Vniax and parametric LFs in Figure [10] 
shows good agreement across all redshift bins for both LF models. 
We find the evolution of the TIR LF is independent of the two LF 
models used (q l = 3.82± ;:j u for model 1 and a L = 3.82± ;?Sj for 
model 2), similar to the results of the 70 and 160 \im LFs. Thus 
the choice of the LF model does not affect the conclusion that the 
TIR LF which shows strong evolution in luminosity. The global 
evolution parameter of the TIR LF is consistent with values deter- 
mined by previous studies, which typically find a L ~ 3.5 ± 0.5 
dCaputi et al]2007l : lLe Floc'h et al]2005l ; lMagnelli et alj2009l) out 
to z ~ 1. In Figure [TO] the TIR LFs are compared with literature 
values in the higher redshift bins. As shown in this figure, our TIR 
LF shown is broadly consistent with the published results in the 
redshift bins [0 - 0.2], [0.2 - 0.4] and [0.4 - 0.8]. We find ex- 
cellent agreement when we compare our 1/Knax L F in the redshift 
bins [0 .3 — 0.451 and [0.6 - 0.8] to the values of iLe Floc'h et al.1 
d2005h and iRodighiero et al.1 d2010h in the same redshift bins. In 
the highest redshift bin, [0.8 — 1.2], our TIR LF shows a good 
agreement with the LF determined bv lLe Floc'h et al.l d2005l) . The 
SWIR E TIR LF h o wever shows di f ferenc es w ith the values of 
ICaputi et al.l d2007h; iMagnelli et alj d2009h and iRodighiero etai] 
Magnelli et al. (2009) investigated the difference between 
their and the Le Floc 'h et al. f d2005h LF values and found that the 
difference appears to be because of the choice of SED library and 
the correlation used to convert L\ to Lir. Thus, this may also be 
the reason for the difference seen in the comparison of this work 
with theirs. Therefore although the 1/Vmax LF values in the highest 
redshift bins are less robust, the parametric TIR LF and its evolu- 
tion conforms to results that have been determined previously. 

In Figure [TTJ we compare our local TIR LF wit h the IRAS re- 
vised b right galaxy sample (BGS) derived local LF of lSanders et al.1 
d2003l) (green filled squares) at z < 0.1, Spitzer 24 [im derived 
LLF of IRodighiero et alj d2010l) (filled blue stars) a t z < 0.3 and 
HerMES SPIRE derived LLF of lVaccari et all fcoich (filled red tri- 
angles) at z < 0.25. The shape of the SWIRE local TIR LF and 
literature values are almost identical while the small difference be- 
tween the local TIR LFs is most likely due to cosmic variance or 
the choice of the SED l ibrary used to calculate Lir. Several authors 
dLe Floc'h etafll2005l : IRodighiero etai]|2010l : IVaccari et al.ll2010h 
have reported best fit parameters of the local TIR LF by fitting to the 
binned LF for the power-law and log-normal LF model of Equation 
[J_4]using a \ 2 minimisation procedure. The slope of the faint end is 
not well constrained and is usually fixed to the local value of ~ 1.2, 
while o values in the range 0.39-0.72 have been reported. The 
logL* (L Q ) range i s ~ 9.24-10.6 and log(/>* (Mpc~ 3 ) range is 
~ -2.00 to -2.06 dLe Floc'h et~ai]|2005l : IRodighiero et alj|20icl: 
IVaccari et al.ll2010l) . The best fit parameters of our LLF for model 
2 are {log L* (L ), a,o, log0* (MpcT 3 )} = {9.87, 1.38, 0.68, 
—2.26}, which are consistent with the published values. For the 
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Figure 10. TIR LF split into redshift bins < z ^ 0.2, 0.2 < z ^ 0.4, 0.4 < z ^ 0.8 and 0.8 < z ^ 1.2. Filled black circles are the 1/V max estimates 
and the black s olid line is the paramet ric Bayesian LF displayed using the best-fit p arameters given in Table [3] at z = 0.1, 0.3, 0.6 and 1. Filled red filled 
circles are fr om Le Floc'h et al. (2005), filled purple upsid e down triangle s are f rom ICaputi et al.l l i2007t) . filled orange filled upside down triangles are from 
iMagnelli etall 120091) and filled blue filled stars are from lRodighiero et all b010l) . The 2 = LF is shown in each panel (black dashed line). 
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Figure 11. The local (0 < z < 0.2) TIR LF determined using 1/V max 
(filled black circles) and the parametric Bayesian method (solid and dashed 
black fine) evaluated at z = 0.1. The filled green squares are obtai ned from 
Sanders et al. (2003) (S03), the filled blue stars are from Rodighiero et al. 
<2010l) (R10~ and the filled red triangles are from IVaccari et all feOld) 
(V10). 



In Figure [ToJ we compare our TIR LF with literature values 
in the higher redshift bins. As shown in this figure, our TIR LF 
shown is broadly consistent with the published results in the red- 
shift bins [0.2 - 0.4] and [0.4 - 0.8]. We find excellent agree- 
ment when we compare our 1/Vmax L F in the redshift bins [0 .3 — 
0.451 and [0.6 - 0.81 to the values of iLe Floc'h et alj fe005t) and 
Rodig hiero et alj d2010h in the same redshift bins. In the highest 
redshift bin considered here, our 1/V max and the parametric LF 
show poor agreement because we are limited by the number of ob- 
jects in these red shift bins. However w hen we compare the LFs wit h 
the values from ICaputi et all d2007l) . and ILe Floc'h etall d2005h . 
we find good agreement between the pa rametric LF and the li tera- 
ture values. Comparing our LF with the Magnell i et al.l d2009l) and 
Rodighiero et alJd201(]|) LF values (orange filled upside down trian- 
gles) shows a slight difference in that t he binned estimates a re lower 
than our parametric and 1/V^nax LF. IMagnelli et al] J2009h i nvesti- 
gated the difference between their and the Le Floc'h et al.l d2005t) 
LF values and found that the difference appears to be because of 
the choice of SED library and the correlation used to convert vh u 
to Lir. Thus, this may also be the reason for the difference seen 
in the comparison of our work with theirs. Therefore although our 
1/V m ax LF values in the highest redshift bins are less robust, our 
parametric TIR LF and its evolution conforms to results that have 
been determined previously. 



double power-law model. [Sanders et all J2003h find best- fit power- 
law indices a = -0.6 ± 0.1 and /3 = -2.2 ± 0.1 and logL* 
~ 10.5 which are at least within the 99.7% errors of the parame- 
ters estimated in our study. 
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;Li R (L ) 0< 2^0.2 



log <£ (Mpc -3 ) 
0.2 < z < 0.4 0.4 < 2 ^ 0.8 



0.8 < 2 < 1.2 



9.2 
9.6 
10.0 
10.4 
10.8 
11.2 
11.6 
12.0 
12.4 
12.8 
13.2 



-1.80 
-2.08 
-2.51 
-2.42 
-2. 
-3.02 
-3.74 
-5.49 

-5.28 



+0.22 
0.49 
+0.18 
0.32 
+0.19 
0.36 
+0.08 
0.10 
+0.07 
-0.09 
+0.07 
0.09 
+0.13 
0.19 
+ 0.G9 
5.49 



+0.31 



,+0.24 
-0.64 
, + 0.11 
-0.15 

-3.09+°-°2 



-2.96Z 
-2.73^ 



-3 
-4.74 



0.00 
+0.08 
-0.10 
+0.17 
0.30 



-3.60 
-3.43 
-3.96 
-4.95 
-5.56 



+0.13 
0.19 
+0.14 
0.21 
+0.08 
0.09 
+0.13 
0.18 
+0.20 
0.38 



-4.53 

-4.57 
-5.42 



+0.24 
0.56 



+0.27 
0.88 
+0.90 
5.42 



Table 5. The 1/V max TIR LF values. 
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Figure 12. Evolution of the space density of normal galaxies (solid and 
dashed blue region), LIRGs (solid and dashed green region) and ULIRGs 
(solid and dashed red region) up to 2 = 1.2. The solid and dashed line 
regions includes 68% of the posterior probability for the double power-law 
and power-law and log-normal LF models respectively. 



6 DISCUSSION 

6.1 Evolution of the Infrared Luminosity Density 

Having determined the global evolution of the TIR LF in the red- 
shift range < z < 1.2, we can now estimate how the comoving 
IR luminosity density (IRLD) and the space density of LIRGs and 
ULIRGs evolves with redshift. This is important because it allows 
us to evaluate the relative importance of IR luminous galaxies and 
their contribution to the CSFRD (associated with obscured star for- 
mation) out to z ~ 1.2. 

We show in Figure [12] the space density of normal galaxies 
(blue regions), LIRGs (green region) and ULIRGs (red region), 
which shows rapid evolution in the number density of LIRGs and 
ULIRGs while the number density of normal galaxies remains 
constant. The number density of LIRGs increases by a factor of 
~ 43^12 for tne double power-law LF model and ~ 37lgg for the 
power-law and log-normal LF model between z = and z = 1, 



which is in agreement with the estimate o f Mag nelli et al] d2009h . 
who find a factor of 40. 

Mag neTli et al.l d2009h estimate an increase in the number den- 
sity of ULIRGs by a factor of ~ 100 at z ~ 1 than in the lo- 
cal Universe, while in Figure [12] ULIRGs increase by a factor of 
~ 316±i3^ for model 1 and for model 2 the factor is ~ 870±^. 
The difference between the values of model 1 and 2 is because of 
the steeper drop-off at the bright end of the LF for the power- law 
and log-normal LF model. Nonetheless, the 68% credible intervals 
shown in Figure[T2]for the two LF models are almost identical ex- 
cept for the ULIRG population, where the error for the power- law 
and log-normal LF model is larger. 

We model the evolution of the space density of ULIRGs as 
oc (1 + z)", with n = 8.30t n'>«, which is consistent with the re- 
sults o f lKim & Sanders! ( fl998h who find n = 7.6 ± 3.2. lGoto et ail 
< l201ll) use FIR data from AKARI to find that ULIRGs evolve with 
n = 10.0 ±0.5, which is sim i lar to the value determined in our 
study. In contrast, Ijacobs et al.l d201ll) have estimated n ~ 6 ± 1 
using a sample of 160 u.m selected sources from Spitzer observa- 
tions of the 1 deg 2 ISO Deep Field region in the LH. Their analysis 
however is based on only 40 galaxies while the current study in- 
cludes 634 sources and therefore the value determined here should 
be more accurate and reliable. 

In Figure [T3l we show the IRLD which is calculated by Air 
= J L<&(L, z\{9}) dL. We only display IRLD determined using 
the double power-law LF model because the 68% credible inter- 
val is almost identical when IRLD is calculated using the power- 
la w and log-norm a l LF model. We follow the method presented 
in lLe Floc'h et al.l d2005T) by determining IRLD for IR luminous 
galaxies (LIRGs and ULIRGs) and "normal" galaxies as a function 
of redshift. The evolution of Qm. is also represented in terms of an 
IR equivalent SFR using the calibration given by Equation[l6] The 
figure shows that our results are in agreement with what has been 
predicted from previou s studies, particularly with the results of 
iRodighiero et alj d2010h . At the present epoch, most (~ 85±°l%) 
of the star formation activity is taking place in normal galaxies 
which have low extinction. As we move to higher redshifts, LIRGs 
evolve rapidly and dominate the star formation activity beyond 
z ~ u.oo_ 05 . 

The model also suggests that, LIRGs and ULIRGs are respon- 
sible for ~ 68±o° per cent of the total IRLD at z = 1.2. At z = 0, 
LIRGs and ULIRGs are responsible for less than 10± |% of the 
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IRLD, while at z = 1, they produce ~ 66lns % of i t. The se val- 
ues are consistent wi th estimates of iLe Floc'h et al.l |2005) while 
Mag nelli et all d2009l) suggest that LIRGs and ULIRGs produce 
less than 2 % of the IRLD . This difference is most likely due to 
the fact that lMagnelli et al.l ( 2009) use the f aint-end power - law in - 
dex, a = -0.6, from the IRAS TIR LF of ISanders etall d2003l) . 
which is steeper than the value used in this study (a = 0.34+QJ2 
for the double power-law LF model) and this would lead to a higher 
contribution of normal galaxies to the IRLD at z — 0. 



6.2 Evolution of the Cosmic Star Formation Rate Density 

In Figure |T4j we show the CSFRD as a function of redshift up to 
2 = 1.2 and compare ou r res ults with the CSFRD estim ates com- 
piled by |Hopkinj ( l2004l) and lHopkins & Beaconj 120061) at wave- 
lengths ranging from the X-ray to radio. We highlight for compari- 
son, in fille d blue squares the SFR de rived from 24 \im observations 
taken from Ruio pakarn et al. (blue filled squares in Figure 

[T4l Overall, our prediction of p S FR show a good agreement with the 
results determined from previous experiments. At z — 0, the CS- 
FRD is 1.04+013 x 10 -2 M yr -1 Mpc" 3 which agrees with pre- 
vious values of the local p SFR that have properly accounted for the 
in ternal extinction of the ga laxy but slightly lower than the estimate 
of Ruiopakarn et al. (2010) wh o find a local p SFR of 1.65 x 10~ 2 
M Q yr" 1 Mpc -3 (assuming a lSalpeterlll95a IMF). T he du st un- 
obscured CSFRD estimated by ISchiminovich et al. (20051) (dot 
dashed orange line) is also displayed in Figure [14] At z — 0, the 



dust unobscured CSFRD is estimated to be ~ 5.01 x 10 M 







yr 1 Mpc 3 and therefore without dust extinction corrections, the 



local CSFRD would be underestimated by ~ 65_3 per cent. Mod- 
elling Psfr as oc (1 + z)® , implies Q = 3. 82+n' 2 R, which is higher 



than the value r eported bvlHopkinsI d2004 (Q = 3.29 ± 0.26) but 
consistent with lRuiopakarn et alj d20ld) (Q = 3.50 ± 0.20). The 
results presented here show that in general, most of the star forma- 
tion over the last 8 billion years has taken place in dust obscured 
galaxies. 



6.3 Future prospects 

In Section [5] we have presented the results of the evolution of the 
far-IR luminosity functions in the SWIRE XMM-LSS and Lock- 
man Hole fields, based on a spectroscopic redshift sample. The 70 
\im and TIR LFs were constructed using the sample selected with 
S70 > 10 mjy and r < 22 and the 160 (im LF was constructed 
using sources with Siso > 60 mjy and r < 22. In this section we 
discuss the limitations of the sample selection and the methodology 
in constructing the FIR LFs and suggest improvements that could 
be carried out in order to obtain better estimates of LFs. 

We chose to study the evolution of the FIR LFs using only 
spectroscopic redshifts in order to reduce the uncertainties related 
to photometric redshifts as stated in Section Q] Some previous 
studies of the evolution of the IR LFs that have utilised photo- 
metric redshift ([perez-Gonzalez et al . 2005; Babbedge ~etalJl2004 
IWang & Rowan-Robinsorj|201(]|) account for the uncertainty in the 
redshift by performing Monte Carlo analysis, where each source's 
photometric redshift is replaced by a redshift drawn randomly from 
a Gaussian distribution centred on the original ph otometric redshift 
and then iteratin g the LF procedure. Few authors dChen et al.l2003l ; 
lAird et alH 2010) have attempted to account for this uncertainty by 
modifying the likelihood function in Equation [10] and including 
the actual redshift probability distribution function for each source 



likelihood (p({Li, Zi}\{9})) rather than a gaussian approximation. 
Therefore including sources with photometric redshifts will not 
only increase the sample size but we will be able to apply a fainter 
magnitude cut in the optical and reach higher redshifts. 

At IR wavelengths, the sample includes sources up to the flux 
limit of the SWIRE survey but with a strict magnitude cut of r < 22 
our sample will not include those IR sources that are also optically 
faint. For the present sample, > 50% of the 70 |J.m sources at 10 
mjy have an optical counterpart that is brighter than 22 mag in the 
r— band and therefore inclusion of the photometric redshifts will 
also reduce the dependence on the spectroscopic incompleteness 
weights Wi. Since future large surveys will discover many millions 
of galaxies, photometric redshifts will play a vital role in any statis- 
tical analysis and therefore the inclusion of this in our studies will 
be one of the next key step. 

The other limitation of the work is that in the parametric 
method, we have assumed a single evolving population. For the 
70 |j.m and the TIR LF the selection functions are calculated as- 
suming an M82 SED template and for the 160 \im LF we have 
used a mixture of M82 and Cirrus SED. As shown in Figure [3] 
we observe a gradual change in the best-fit SED type as we move 
out to higher redshift, with low redshift sources fit with a cirrus 
dominated SED and the higher redshift sources fit with a starburst 
dominated SED. This implies that the starburst component of the 
IR SED must evolve much faster than the cirrus component and 
therefore this evolution must be taken into account when studying 
the evolution of LFs using parametric methods. Therefore studying 
the evolution of each SED component will be the second step in 
studying the evolution of the FIR LFs. 



7 SUMMARY AND CONCLUSION 

We have presented a new observational determination of the FIR 
LFs using 70 vim selected sources from the SWIRE survey with 
spectroscopic redshifts. The primar y sample was selected f r om th e 
photometric redshift catalogueue of Ro wan-Robinson et al. I d2008l) . 
which contains over 1 million IR sources estimated by combining 
optical and IRAC 3.6 and 4.5 urn photometry. We computed the 
rest-frame 70 [im and TIR LF using sources with SVo > lOmJy and 
r < 22 and studied their evolution out to z = 1.2. The evolution of 
the 160 \im LF was determined using Sieo > 60mJy and r < 22 
out to z = 0.5. 

We use the multiwavelength optical to IR data to model the 
SED for each source in our sample to estimate the rest-frame 
monochromatic luminosities to determine the 70 and 160 \im LFs 
using Bayesian parametric and the 1/Vma.x methods. The work pre- 
sented here is an improvement on earlier works because we use 70 
\im data rather than rely on extrapolations from 24 [im. In addition, 
the Bayesian method is a new approach to FIR astronomy, which 
can be used to further the study of IR LFs. For example we can 
include photometric redshift data and properly account for the in- 
dividual photometric redshift uncertainty of each source to investi- 
gate the evolution of the FIR LFs. The study of the evolution of FIR 
LFs using photometric redshifts will be presented in a future paper. 
In our analysis we have corrected for the optical and IR selection 
biases and incompleteness of the spectroscopic sample, which have 
allowed us to accurately construct the LFs. 

The parametric Bayesian FIR LFs were determined using two 
LF models; a double power-law and a power-law and log-normal 
models. Comparison of the two models with the binned 1/V ma x 
estimates showed that both LF models provide a good fit to the 
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Figure 13. Evolution of the IRLD up to z = 1.2 (solid black line region) and the contribution by normal galaxies (solid blue line region), LIRGs (solid 
green line region) and ULIRGs (solid red line region) determin ed usin g the double power-law LF model. The black solid region includes 68% of the posterior 
probabi lity. The red lin e region is taken from iLe Floc'h etafl l2005ri . orange line region is from iMagnelli et al.1 d2009h and blue filled stars are taken from 
iRodighiero et alJfeOlOl) . 
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Figure 14. Evolution of the CS FRD up to z = 1.2 (solid black line region). The solid black line region includes 68% of the posterior probabil ity. The grey 
points are taken from[Hopkins (2004) representing the CSFRD from various estimat ors. The red filled squares are CSFRD estimates from Ruiopakar n et al.l 
120101) . The orange dash-dot line is the UV dust extinction uncorrected CSFRD from Schiminovich et al. ( 2005 ) and the light blue dashed line represents the 
total CSFR D defined as the sum of the IR and UV dust extinction uncorrected SFR densities. The shaded grey region is the 3<r best-fitting regions CSFRD 
from lHopkins & Beaconj 120061) . 
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70 Jim and TIR LFs. For the 160 vim LF however, the power-law 
and log-normal LF was found to underestimate the number density 
of sources at Li6o > 10 1 Lq when compared to the binned LF 
estimates, implying that the double power-law LF model provides 
a better description of FIR LFs. 

The evolution of FIR LFs was modelled using a pure luminos- 
ity evolution model. The rest-frame 70 vim LF was found to evolve 
rapidly in luminosity with a L — 3.39^q'22 for model 1 and a L 
= 3.411q25 for model 2. Thus the choice of LF model does not 
have a large effect on the evolution of the 70 vim LF. The 160 vim 
LF was found to evolve with a L = 5.73^q'45 for model 1 and q l 
= 5.621q'45 for model 2. The fa ster rate of evolution o f the 160 vim 
LF is consistent with the study o flTakeuchi et~al] ( l2006h which finds 
the ISO 170 vim LF evolves with a L = 5.OI05 . Since cooler galax- 
ies are detected more effectively at 160 vim than at shorter wave- 
lengths, the strong evolution implies the presence of large amounts 
of cool dust at higher redshifts. 

Finally, we used the best-fit SEDs to estimate total IR lumi- 
nosities of each source to derive the TIR LF to z = 1.2. The 
SWIRE TIR LF showed good agreement when compared with lit- 
erature values estimated from Spitzer and Herschel data across all 
redshifts bins. The TIR LF evolves with a L = 3.821q'2q for model 
1 and q l = 3.82±° for model 2 corroborating the conclusions 
of previous results from Spitzer 24 vim studies which find strong 
luminosity evolution. The TIR LF was then integrated to calculate 
the IRLD out to z = 1.2, which confirms the rapid evolution in 
number density of LIRGs which contribute ~ 66I5 per cent to 
the IRLD, and hence the CSFRD at z = 1. The results presented in 
this chapter, based on 70 vim data confirm that the bulk of the star 
formation at z — 1 takes place in dust obscured objects. 
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